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1 Introduction 

> : 

This paper is concerned with the linear theory of the generalized thermodynamics 
e$ • for the micropolar porous media. 

The micropolar elasticity theory has been established by Eringen [1-4] and is of 
fundamental interest in a lot of cases for which the classical theory of elasticity is 
inadequate. Other related results have been obtained in [5-7]. 

A general theory of materials with voids has been developed by Cowin and Nun- 
ziato in [8-9], by using the concept of distributed body introduced in |I0| . In this 
theory, the matrix material is elastic and the interstices are void of material. In this 
context, the bulk density is written as the product of two fields: the material density 
field and the volume fraction field. The intended applications of this theory are to 
geological materials (like rock and soils) and to manufactured porous materials (like 
ceramics and pressed powders). 
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In a previous paper | TT| , we studied the thermodynamic theory of micropolar 
porous materials by using the theories developed in [1-4], || and in [8-9]. This study 
was based on a generalized theory of thermoelasticity including the heat flux vector 



in the set of the constitutive variables. In [11 1, we also derived the corresponding 
linear theory. 

In the present work, starting from the basic equations established in |TTJ, we state 
two variational characterization (of Biot and Hamilton types, see [12], [13]) of the 
initial-boundary value problem. In final part, we also derive a reciprocal relation 
for the problem in concern by using convolution in time (fundamental results, in 



this context, have been obtained by D. Grafli []I4" |) and Laplace trasform. 



2 Basic Equations 

Throughout this paper, by Q and dQ we shall denote the region and its boundary 
of the physical 3-dimensional space (= 3? 3 ) occupied by an elastic solid in a given 
reference (unstressed) state and / = [0, +oo). We refer the motion of the body 
to a fixed orthonormal frame in 3? 3 . We shall denote the tensor components of 
order p > 1 by Latin subscripts, ranging over {1,2,3}. Summation over repeated 
subscripts is implied. Superposed dots or subscripts preceded by a comma will mean 
partial derivative with respect to the time or the corresponding coordinates. 

On the basis of the linear theory established in |Tl| , the behaviour of a micropolar 



porous thermoelastic body is governed by the following local balances of momentum, 
equilibrated force and energy 

Tjij + Pfi = phi, 

(1) 

hi,i + g + pi = pX (p, 
p(6 T] - r) = q ih 
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where 

Tij, Mij stress and couple stress tensors; 

hi, g equilibrated stress vector, intrinsic equilibrated body force; 
T], qi entropy and heat flux vector. 

These are the dependent variables of the theory, needing a constitutive equation. 

Further, p is the bulk mass density, the body force, £{ the body couple, J^- 
the micro inertia tensor, X the equilibrated inertia, I the extrinsic equilibrated body 
force, r the strength of the internal heat source and 8q the (constant > 0) tempera- 
ture in the reference state; €i rs is the alternating symbol. 

Finally, Ui, (pi and <p, along with the temperature 9 (measured from 9q), are the 
thermokinetic variables of the theory: «j represents the displacement, fa the micro- 
rotation, (p the change in volume fraction from the reference configuration. 

The constitutive equations for ([[]) are |JTTJ 



(2) 



= CijkiEki + Gijki&ik + Hijip + Hijktp^k + Aij8, 

= GkUjE k i + r^ki^ik + Pij(p + Pijk^p,k + Gij6, 

9 = -HijEij - Pij&ji -atp- ditfj - b9, 

hi = H jki Ej k + Pjki&kj + aw + Dijifj + 7j0, 

pr] = -AijEij - Gij&ji -bip- ji(p }i + c9, 

(l + r|) ffi = KyOj. 
Here, J5y and are the deformation and wryness tensors defined by 

Eij = u jt i - (,ji.Oi,. fyj = Oi. r (3) 
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We shall consider a homogeneous body; in this case, the coefficients in (|2|) are all 
constants; it can also be proved (Tl| that the following symmetry relations follow 



and that the conductivity tensor is inversable. 

Note that applying the operator [\ + rjQ to the energy equation, by we 
get an evolution equation for 9 of hyperbolic type 

+ r^j [pr - p9 f)] = 0. (5) 

By insertion of (0, into (|TJi,2,3, H), we obtain the complete system of (coupled) 
field equations for u i: (pi, <p and 9. 

Let us denote four pairs of disjoint and complementary subsets of the (smooth) 
boundary dQ by Sj+i} with i — 1, 3, 5, 7 and denote the outward unit normal 
to dfl by rij. 

To the system of field equations (^1,2,3 , @ , 0) we add the following initial-boundary 
conditions 

Ui = iii = 4>i = 4>i = ip = if = 6 = qi = in Q x {0}, (6) 



Ui = u* on Y>i x /, 
k = & on S 3 x /, 



TjiUj = t* on S 2 x /, 



MjiUj = m* on S 4 x /, 



Lf = Lf* OI1E5 X /, 



hjUj = h* on Eg x /, 



(7) 



9 = 9* on E 7 x /, 



qjUj = q* on S 8 x /, 



4 



where right-hand terms stand for (sufficiently smooth) assigned fields; along with 
fi, £i, i and r, these are the data of the mixed problem considered. An array field 
(ui, (pi, <f, 6) meeting all equations (0i,2,3, @, H 0, 0), for some assignment of the 
data, will be referred to as a (regular) solution of this. 

The uniqueness of solution to the above problem has been established in [ITI 



3 Variational principle I 



In this section, following we establish a variational theorem of Hamilton type. 



We introduce (see also |TTJ) the kinetic energy K and the free energy ip 
pK = \p(uiUi + J ij 4>i4> j + X (f 2 ) 

pip = \CijkiEijEti + ^rijM&ji&ik + \a(p 2 + ^Dijtpjipj - \cd 2 + 

(8) 

+ \Bijqiqj + GijkiEij&ik + HijEijip + Hij k Eijip tk + AijEij9+ 
■ I'.^j.r + Pijk&ji¥,k + Gifted + a^y^ + hp9 + %(p }i 9. 
We can prove the following 

Variational principle I Let t± and £2 be arbitrary instants of time in I and 
assume that 

i. the symmetry relations (fjj hold; 

ii. Kij and Jij are symmetric tensors. 

Then, for arbitrary variations in Ui, (pi, ip, 9 such that 

5ui = 0, 5(pi = 0, 5cp = 0, 59 = 0, in Q x {t x , t 2 }, (9) 
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5ui = on Ei x /, 
5 <p = on S 5 x J, 



5 0i = on S3 x /, 
5# = on£ 7 x/, 



(10) 



and the functions fi, Ei, I, rj, r, t* , m* , h* , q* , t\ and t2 being kept unchanged, the 
following variational equations 

5 f 2 dt f p(K-ij-r]T+ ^-B ijqiqj + f iUi + ^ + tip) dV + 
h x Jn v Ip J 

(11) 

+ / t*Ui da + frii4>i da + h*ip da =0 

-/So JY>a -^Sr 



5 / dt 

Jti 



-KijTtTj dV-(l + rp)\ p( V T - r)T dV + / q*T da 



(12) 



are satisfied if and only if are verified the field equations (^23) an d (0/ ^ n the 

d 

above equations we put T = 9 + 9 (the absolute temperature) and p = — . 

at 

Proof. We can remark that, under the conditions (§), the variation of the ki- 
netic energy is 

5 J 2 dt J pK dV = - J 2 dt J p(uidui + J ij 4> j d(f)i+ X <p6<p) dV. (13) 



Moreover, using (^2, @, 0, f|, 0; @> 0) an d the divergence theorem, we can obtain 
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^ It dt In (f^ + ^ ~ \ Bijqiq I) dV = 

= f 2 dt f {TjiSEji + MjiSWij - gStp + hiStp^dV = 

= j dt f {-TjijSui - (Mjij + e irs T rs )54>i - (h iti + g)8tp} dV + 



(14) 



+ / t*5u{ da + m*5(f)i da + h*8ip da 



Taking into account (|1^) and ([13]), the variational equation ( |Tl~D can be rewritten as 
follows 



ta 

dt 
ti Jn 



(Tjij + pfi - pui)5ui + {M^j + e irs T rs + pl { - pj^^b^ 



+(hi,i + g + pi - pX (p)Scp 



dV = 



This equation holds if and only if (|l|i,2,3) are satisfied. 
From (BO and (M) it follows that 



S I dt -K^TiTj dV = I dt 



f 2 dt\(l + rp) f q*59 da- f Kij6 fi j86 dV 



where p has been treated as a constant [0 . 



Thus, fll2|) can be written as 



t dt f {KijO^ + p(l + rp){r - Ti])} 56 dV = - f {r]T5T)\\\ dV 
Jti J n Jn 



After the linearization, from (B|) we conclude that 



f 2 dt f {KijO^ + p(l + rp){r -9 r))}56dV = 
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This equation holds if and only if @ is satisfied. 



4 Variational pinciple II 

Following ||12|| , we introduce the entropy flow vector Sj by 

Si = qi 

and by (|1|4, g;) we get 

s M = -AijEij - Gij&ji - bip - 7i0,j + c6 - ^-r. (15) 

fro 

In the following, we shall consider prescribed f\, i and r in Q, and Ui on E 1; 
on E2, 4>i on S3, Mij on S4, ip on E5, hi on Eg, # on E7 and on Eg. Moreover, we 
also consider Sj assigned on E 8 . 

If we assume that is symmetric and denote its inverse by Ky (= Kjij , we set 



D = ^Jjv{p + Tp 2 )K ijSlSj dV 



so the variation of this functional due to variation of is 



6D = J q (p + Tp 2 )K ijSl 5 Sj dV. (16) 



Here, p is used as in section 3. 

We now introduce the following functional 



V = ( (pip + pr]0 - -B^q^) dV. 
If we consider variations of V due to variations of Ui, (pi, <p, 9 and use (§2), the 
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constitutive equations and the symmetry relations (^), we obtain 

sv = / {cese + (x tJ - .\, ! o)6i-:, J + (m^- - g^)c%+ 

(17) 

-((/ + W)^+(A i - 7i 0)^, i }dV. 
Moreover, we define the following functional 

G = — t*Uido — I m*(pida— / h*ipda— / 0*SiHida 

If 5(f>i , 5ip, 5si vanish on E 1; E 3 , S 5 , E 8 , respectively, by (|7|) we can write 
5G = — J (Tji5ui + Mji5(pi + hj5(p + 85sj^jrij da; 

and, from the divergence theorem and (|3|, f|, [T5|), we get 

5G= - I {(T tJ - .\, J 0)6l-:, J + (Mij - G tj e)6^ + (h t - 7i 0)<% + c656+ 



(18) 

—b68(p + Tji jdui + (Mjij + e irs T rs )8(j)i + hjj8(p + 9,j8sj} dV. 
Finally, we put 

F = L p { {l p2ui ~ fi ) Ui + {l p2j ^ j ~ £t ) & + G x p2 ^ ~ £ ) V ] dv > 

so, assuming the tensor Jij symmetric, the variation of this functional is 

5F = J p{(p 2 Ui - fi)Sui + {p 2 Jij(j)j - ti)Hi + (X P 2 <p - £)5<p} dV. (19) 

Thus, we can obtain the following variational theorem (of the Biot type) 



Variational pinciple II Let us prescribe fi, i; L , £ and r in £7, and T^j, My, hi, 
and 9 on E2, E4, E6, E7, respectively, and assume that 

i. the symmetry relations hold; 

ii. Kij and are symmetric tensors. 

Then, for every 5ui (= on Ex), 5(pi (= on E 3 ), 5<p (= on E 5 ), 5si (= on 
E 8 ), the following equation 

5H( Ui , (j)^ <p, 9) = 5V + 5G + 5F + 5D = 

yields the field equations (^1,2,3 j HI)- 

Proof. By (0, |T7|, [IB], [H5]) we can write 

SH(ui, (pi, if, 9) = f n {(pui - Pfi ~ Tji d )b~Ui+ 

+ (pJij<Pi - ph - M jid - e irs T rs )5(pi+ 

+(pX(p-p£-g- h it i)5ip+ 

+ [0 O (1 + rpiKijSi - e ;j ]5 Si } dV. 

Vanishing of this expression for arbitrary 5ui, 5 (pi, 5tp, 5si of course implies the 
field equations [l]i,2,3; further from the last term in the integral above, we get 

O (1 + Tp)KijSi = 9 j. 

If we multiply the above equation by Khi and derive with respect to the i th coor- 
dinate, we can easily deduce the field equation (g) from (^). This completes the 
proof. 
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5 Reciprocity 



Let us denote by a * b the time convolution of the scalar fields on f2 x / 

(ai * a 2 )(x, t) = / ai(x, £ — r)a 2 (x, r) o?r. 
Jo 

Let us consider the body subjected to two different sets of data 

T («) = |jW ^ ^ ^ r (a) ? u * (a) ^ t * (a) ? ^ (a) ^ ^ (a) ? ^* (a) ? ^* ( a ) ? ( a ) ^ ? , (a) j 

and null initial conditions and let be 

M (a) = («i a) ,#,V (a) i e (a) ) a = 1,2, 

the corresponding solutions. Moreover, we define h[ a \ g^ a \ i]( a \ qf^ by 

means of (0) for each a = 1,2. 
We prove the following theorem 

Reciprocal Theorem Let U^ a ' be solutions corresponding to different sets of 
data T^ a > (a = 1,2), and assume that 

i. the symmetry relations (fjj hold; 

ii. Kij and Jij are symmetric tensors. 
Then, the following relation holds 

lap = 1p a (20) 

where 
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T, 



a(3= IP 

n 



^ ii* {a) * T^n 3 da + / t* {a) * da+ 



i* {Q) * Mjf ^ da + i m* {a) * </>f da+ 



s 4 



^ 0* {a) * hfrij da+ h* (a) * ^ da 



+ J_ / 0* («) * g f ) n . da- - [ q* (a) * # (/3) da 
6*0 js 7 6>o Je 8 



for a, /? = 1, 2. 



Proof. We put 



T 

** ( 



(\ 1 



where a is Laplace transform with respect to t of the function a: d(s) = / °° e~ st a(t)dt. 
If we take Laplace transform of ([3], 0) under the hypotheses i, ii and the initial con- 
ditions (10), we deduce that 



12 



7~£af3 — 7~£f3a 



-9 s \A l3 (6^Ef - 6We\?) + G^O^P - 6^ y )+ 



so 



(21) 



If we integrate this equation over Q and use the divergence theorem and boundary 
conditions, we get 



ps 



r«f ) +r# ) +^ ) ^ ) -F ) «r ) + 



jWffi _ ^)^(a)l _ n r («)0(P) _ r (P)0(a)l d y + 



P_ 

9o 



Si 



^2 



+ / S (Mjr } 0l (/3) - Mf4>: [a) )n 3 da+ [ s{m^W - m*^) da+ 
+ f s{hf(p*W - hfip*^)nj da+ f sCh* ia) v m - h*^(p^) da+ 



Eg 



E 6 



9q Jt, 7 &o J^b 



(22) 



Finally, we obtain the desired equation (gOj) by inverting (P2|) and using the convo- 
lution theorem for Laplace transforms. 
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